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D , Abstract. The equation of state is investigated for an Ising-like model in the 

framework of collective variables method. The peculiar feature of the theory is that 
a non-classical van der Waals loop is obtained. The results are compared with ones 
of a trigonometric parametric model, in terms of normalized magnetization, M, and 
field, H. 
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1. Introduction 

Recently, in [1] it was shown that considering the system of Ising spins on a simple 
cubic lattice by the collective variables (CV) method (H [3], a non-classical van der 
Waals (vdW) loop can be obtained. In the present paper, the shape of the loop is 
investigated. A comparison is made with similar results of [I] , where the van der Waals 
loop was obtained for a system with a scalar order parameter using a trigonometric 
parametric model for scaling behaviour near criticality. The discussion on whether such 
a loop exists naturally and has or not any physical manifestation can be found in [5]. 
Our purpose in this paper is to compare results from CV theory with others. 

We consider a system of N Ising spins on a simple cubic lattice of spacing c. The 
Hamiltonian of such a system is well known 

H = ~l'£*(r i j)a i a j -H y E l a i . (1) 

i,j i 

Here the spin variables o~i take on ±1, H is the external field, and Qfrij) is a short-range 
interaction potential between spins located at the i-th and j-th sites of separation r^. 
The interaction potential can be chosen in the form of exponentially decreasing function, 
= const ■ exp (— r„/6), with b being an effective range. 
The partition function 

Z = 5>-'* (2) 
M 

where (3 = (k^T) -1 is the inverse temperature, can be written in terms of collective 
variables pk (2J [3]. In "p 4 - model" approximation, the explicit form for such a 
representation is as follows 



J (dp)^ exp a lV /iVop d {k)f>*P-* ~ 



^J^O 1 Y Pki • • ■ Pk 4 <5ki+...+k 4 



keBo 

(3) 



Here the quantity d(k) contains the Fourier transform of the interaction potential 

d(k) = a 2 + /3$ (0) $-/?$( k). (4) 
For $(&) we use the so-called parabolic approximation 

[% = $(0)$ for B Q < k < B. 

Strictly speaking, the wave vector k takes on values from the first Brillouin zone 

7T 27T Tl ■ 

B = {k = (k x ,k y ,k z )\ki = 1 -rf;^ = 1, 2, . . . , N t ] % = x,y,z}. (6) 

c c Ni 

In the following, however, we will keep to the spherical approximation for the Brillouine 
zone so that B = tt/c is the boundary of this zone, and B = ir/(cs ) = n/c is the 
boundary of B . The discussion on the choice for s and $ can be found in [B]. In 
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general, Sq should depend on the ratio of the effective interaction range b to the lattice 
constant c. However, in the present calculation, we fix sq = 2 and b/c = 0.3. This gives 
$ = 0.27492. The quantities a* from ([3]) are expressed as 

ai = s d J 2 h, a 2 = 1 - Sq d , a 4 = 2sQ rf , (7) 

where d — 3 is the space dimension, and h = /3"H is dimensionless field. In ([3]) the 
collective variables pk with Bq < k < B have already been integrated out so that 
Nq = Nsq 3 is a number of variables remaining to be integrated. 

Calculation of the partition function ((3]) is performed according to Yukhnovskii's 
method [2] . It is based on the idea of step-by-step integration of the partition function, 
over the subsets of collective variables, first with B\ < k < Bq, then with B 2 < k < B\, 
and so on, while averaging the Fourier image of the interaction potential on each step. 
Here B n = B n _i/s = Bq/s 11 , where s is the renormalization group (RG) parameter. 
This is equivalent to the Kadanoff scheme of constructing spin blocks [7J E\- Every 
time when integrating over a subset of CV, a factor - denote it Q n - appears in the 
partition function. On performing step-by-step integration of the partition function over 



rip subsets, one arrives at 



Z = Z [Q(d)] No ( UQ n )Z hGR . (8) 



n=l 



The partial partition functions Q n of the n-th level are characterized by a set of 
coefficients a^, d n (0), for which the recurrence relations (RR) hold The RR 
have a fixed point as a partial solution. Because of that, the quantity n p is chosen from 
the requirement that for n < n p the RR can be linearized near the fixed point. In this 
case, the system possesses the RG symmetry and is said to be in the critical regime 
of the order parameter fluctuations. The quantity n p is called the exit point from the 
critical regime. 

The analytic expressions for the exit point were found for the limiting cases, h = 
and r = : 

In 1^1 

n p = — : — 1, at r = 0, 



In Si 
In If I 



at h = 0, (9) 



lnE 2 

where r = (T — T c )/T c is the reduced temperature. 

Here E\ = 24.551 and E 2 = 8.306 are the eigenvalues of the matrix of the RG 
transformation linearized near the fixed point of RR, 

d/2 

f=^T, h= S -^-h, (10) 

Jo h 

where fo = 0.5 defines the fixed-point coordinates, cik = 1.155, h = 0.760. In general 
case, the expression for n p cannot be obtained analytically. For example, in [10] this 
quantity was computed numerically as a solution to a certain equation. In any case, n p 
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should satisfy in the mentioned limiting cases. Based on this requirement, in [TTJ 
the expression was constructed for the exit point in the form 

(±) = _ lnffi + )»W 2 ) _ 

' 2 In Ei ' 11 

where some temperature fields, = It] 135 and h*~^ = Iti] 135 , are introduced, (3 and S are 
the critical exponents of magnetization. In our approach, (3 = 0.3024 and 5 = d + 2 = 5. 
The signs "+" and "-" are related to T > T c and T < T c , respectively. In the following, 
we will mainly omit superscript ±. Finally, T\ = —E 2 °t, where Uq denotes the difference 
between and n^T'. In [12] n = 0.5 is chosen to recover the universal ratio of critical 
amplitudes for the correlation length, = 1.896(10) [T3] . 

In (jSJ), the quantity ^lgr is still remaining undefined. Here, we shall not delve into 
considering this quantity. A detailed explanation of how to compute it can be found in 
[j~2| [14"] . We just recall that it is expressed in the form 

Z LGR = Z G Z TR -e NE °^\ (12) 

where the contribution e NE °^ is the most important because it is due to the collective 
variable po- 111 the present research, main attention is paid on this part of the partition 
function. The quantity Eq(o~) has the form 

E = e hCh 2 + h^ 2 )^j _ e2 (/,2 + /,(±) 2 )^. (13) 
The following notation is used 

e = a s~ 1/2 , 

e 2 = 2 a o s ~ 3 ( r « P +2 + j^sla 2 u np+2 ), (14) 

where 

r„ p+2 =/3$(0)/ (-l±^c), 

u np+2 = (/3$(0))Vo(l ± $E 2 H C ), (15) 

with H c = f(h 2 + /j,(±) 2 ) -1 / 2po . The quantity <To is the solution of the following cubic 
equation 

a 3 + po-o + q = 0, (16) 
with the coefficients 

_o T n _|_2 
p = Q Sq — 2— 

Un p +2 

q = - Qsf /2 .^ 2 ^ (17) 

«n p +2 (fc + h (±) y/2 
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2. Equation of state 



Based on (JBJ) and (112I) . the Gibbs free energy can be expressed as a sum of three terms 

F(r, h) = -k B T \nZ = F & + + F (±) . (18) 

Here, the term F a is the analytical part of the free energy and does not affect the critical 
behaviour. The term F^ is expressed as 

F s W(r, h) = -A; B TiV 7 ( ± )(/ i 2 + h^ 2 )^ , (19) 

where jf 1 includes contributions from the critical regime and from the limiting Gaussian 
regime (inverse Gaussian regime in the case of T < T c ). The explicit expression for it 
can be found in (5.6) of [13]. Finally, the quantity Fq^ from (|T8|) is 

Ft\r,h) = -k B TNE (a). (20) 

This contribution to the Gibbs free energy is due to the collective variable po, which, as 
is known from the theory of collective variables [2], is related to the order parameter. 
Therefore, F^ is the free energy of ordering. The order parameter of the considered 
system - the magnetization - is then calculated by means of the thermodynamic formula 

M„ = = (21) 

N dn dh K J 

that leads to the expression 



3/2 

e o( 1 + ~Z -Z-772 ) - 7 e 2-,_ - 2 



, ■■» / 1 h 2 \ 6 s n h 

M = (h 2 + h 2 c )^+5) ■ 1 1 



de2 ' {h 2 + h^yi 



5/, 2 + /l (±) J 5 ho (h2 + h &) )l/2 

(22) 



dh 

The magnetization f l22|) can be rewritten in more compact form 

M = a 00 (h 2 + /4 ±)2 )^ , (23) 

where 5 = 5 is the critical exponent describing the field dependence of magnetization 
and 



o"oo 



1 h 2 \ 6 s 3 /2 h 

— — Q ) — ~&2~, Z 2~ 

5^2 1 a(±) / 5 hn (u2 1 *,(±) ' 



^2 + / a ft (^2 + )i/ 2 



de 2 \ ,r, . , ^2 X 



1/2 



In figured] the magnetization (1221) is presented as a function of the external field. As 
is seen from the picture, the proposed approach gives a van der Waals loop, however with 
non-classical exponents. In order to compare our results with ones obtained by different 
methods, we appeal to [3], where a non-classical van der Waals loop was obtained 
with the help of a trigonometric parametric model for scaling behaviour near criticality. 
For this reason, in the next section the main features of the conventional parametric 
presentations of the equation of state are presented as well as the novelty suggested 
in [4]. 
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Figure 1. Magnetization (|22|) as a function of field below the critical temperature, 
t = -0.001. 

3. Parametric representation of the equation of state 

The way of representation of the critical behaviour with the aid of new parameters, r 
and 9, was firstly suggested in [15]. The meaning of the parameters is that the "radial" 
variable r defines the distance from the critical point on the "temperature- field" plane, 
and the "angular" variable 9 defines the path along a contour of constant r. Today, the 
more traditional form of the parametric representation is as follows 

M = r* 3 m(9), h = r? 5 l(9), r = rk(9), (25) 

where, in terms of magnet, M is the magnetization, h is the magnetic field, and r is the 
reduced temperature. 

The function k(9) is usually chosen in the form k(9) = 1 — b 2 9 2 , which means that 
9 C = ±1/6 corresponds to the critical isotherm, for h > and h < 0, respectively. The 
function 1(9) is chosen in the form of a polynomial with only odd powers of 9. Zeros 
of 1(9) correspond to the absence of the field, h = 0. Particularly, 9 = specifies the 
critical isochore, h = and r = r > 0. Next zeros of 1(9), 9 = ±9%, correspond to 
two sides of the coexistence curve. 1(9) is usually written in the form 1(9) = loh(9) 
so that dl(9)/d9\ e=0 = l'(0) = l , and Z is a normalization constant. And, finally, 
m(9) = mo(9 + c9 3 ). If c = 0, the transformation ( 125]) is referred to as the "linear" 
model, and in the case of c ^ it is called the "cubic" model. 

As can be easily seen, such a representation fails in the two phase region. For this 
reason, in [5], H] were suggested new, the so-called trigonometric parametric models that 
would naturally extend to the two phase region and recover a van der Waals loop. Here 
we present the main features of the extended sine model [1] - the more preferable of the 
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trigonometric models. 

First of all, in this model a transformation of type (125]) touches the temperature, 
magnetization and the thermodynamic potential A(M, T): 

r = rk(6), M = r^m(6), A = r 2 - a n(6), (26) 

where the angular functions are chosen in the form 

k(9) = 1 -2b 2 [l- cos (q6)]/q 2 , 
m{6) = m sin (q9) / q, 

n{6)= Z o m o [l+£ c ^(0)]- (27) 
i=i 

Here the coefficients b, q, and q are regarded as adjustable parameters. In such 
circumstances, for the field angular function 1(9) one has 

im _ j k(9)n[(9)-(2-a)k'(9)n l (9) 

[> k(9)m' l (9)-^'(9)m 1 (9) ' 1 ' 

where we have introduced the following notation 

ni{9)= ^ = 1 + E C » 

mi (0) = H^l = sin ( q 0)/q. (29) 
m 

The van der Waals loop is obtained in terms of normalized versions of M and h, 
namely 

M = M/M (T), H = h/[M (T)/xt(T)\, (30) 

where Mq(T) is the "spontaneous" value of the order parameter, and Xo{T) is the 
zero-field susceptibility above T c , Xo(T) w C + t~ 1 . 

In figure El the dashed line is the van der Waals loop for the following values of the 
adjustable parameters (see formula (6.2) in [1]): 

6 = 3.81857, q = 3.03134, ci = 1.71886 

c 2 = -5.97266, c 3 = -0.207431, c 4 = 0.0882037. (31) 

Note, that for such a choice l'(0) ^ Iq, in contrast to either the "linear" or "cubic" 
model. 

In the same figure [21 the solid line is our results based on (122]) . Regarding the 
approximations, "p 4 -model", and zero value of the small critical exponent rj = 0, we see 
good qualitative agreement in results, especially in the regions of h and of h ^> h c , 
H > 2. Some qualitative discrepancy is observed for the intermediate values of H, i.e., 
for h ~ h c . In this domain, with increasing H our curve initially withdraws from the 
dashed curve, further approaches and intersects it and then moves along with the dashed 
curve, below it all the time. Such discrepancy is connected with the choice of functional 
form of n p (fTTj) . Although it gives correct values in the limiting cases, see Qj, however 
it seems to fail in the intermediate region, h w h c . Therefore, to improve our results we 
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Figure 2. Scaled van der Waals loop and magnetization above the critical 
temperature. The dashed curves represent results obtained within the extended 
sine model. The full curves represent the results of the collective variables method, 
equation (j2"2")l . 



need a somewhat different functional form in this value range. Another way to do so is 
to compute n p numerically, but we do not lose hope to solve the problem analytically 
and will attempt to find more appropriate expression for the exit point in a future work. 

With regard to the loop, we observe the loop is wider in CV theory. The 
spinodal value of the magnetization in CV, M = 0.55 is less then the corresponding 
value M = 0.697 from j5], which means that the CV theory, at least in "p 4 -model" 
approximation, gives wider metastable domain. 

4. Conclusions 

In conclusion, this work is the first attempt to investigate the van der Waals loop 
using the collective variables method, which is essentially microscopic approach. This 
investigation is important because there is lack of non-classical theories that give the 
van der Waals loop (see, e.g., [16], Sec. 3.4.4). Some quantitative disagreement of 
our result in comparison to ones obtained in [I] can be associated with the model 
approximation, since present calculations are carried out using the simplest non- 
Gaussian approximation, i.e., //-model. The obtained result can be improved both 
in a formal way, by appropriately choosing the functional form of the exit point n p , and 
in a conceptual way, to which can be attributed (a) investigating the system in a higher 
approximation, i.e., p n -model with n > 4 (for n = 6 see [17J), (b) including corrections 
to scaling [18], and/or (c) accounting for averaging the interaction potential [19], the 
latter resulting in rj ^ 0. All of these will be the scope of a forthcoming paper. 
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